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ABSTRACT 



Metal-rich (red) globular clusters in massive galaxies are, on average, smaller 
than metal-poor (blue) globular clusters. One of the possible explanations for 
this phenomenon is that the two populations of clusters have different spatial 
distributions. We test this idea by comparing clusters observed in unusually deep, 
high signal-to-noise images of M87 with a simulated globular cluster population 
in which the red and blue clusters have different spatial distributions, matching 
the observations. We compare the overall distribution of cluster effective radii 
as well as the relationship between effective radius and galactocentric distance 
for both the observed and simulated red and blue subpopulations. We find that 
the different spatial distributions does not produce a significant size difference 
between the red and blue subpopulations as a whole, or at a given galactocentric 
distance. These results suggest that the size difference between red and blue 
globular clusters is likely due to differences during formation or later evolution. 



Subject headings: galaxies: individual (M87) - galaxies: kinematics and dynamics 
- globular clusters: general 



1. Introduction 



Globular cluster p opulations within many types of galaxies have colour distributions 
which are bimodal (e.g. IZepf fc Ashman! Il993t iLarsen et al.ll200ll : iHarrisI l2009bl : iPeng et al. 
20061 ). This split in colour is assumed to be due to a wide range in cluster meta l licity. 



with red clusters being metal rich compared to blue clusters (e.g. iZepf &: Ashman! Il993 
Brodie fc Straderll2006[ ). The fact that this bimodality is seen in all but the smallest galaxies 
indicates that globular clusters may share a common formation mechanism and have similar 
stages of evolution from one galaxy to the next. 
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An interesting difference between these two sub-populations has emerged from a long 
list of observational studies which find that blue globular clusters have effective radii {r^) 
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There exist several competing ideas regarding the source of this size difference. iLarsen fc Brodie 



( I2OO3I ) suggested that the size different is not intrinsic, and is due to a projection effect of the 
different spatial distributions of the red and blue subsystems. Many studies show that the 
red sub-population t ends to be more centrally concentrated in the host ga laxy than the blue 



sub-population (e.g. iLarsen et al.ll200ll : iForbes et al.ll2006t lHarrisll2009al Jbl) . From basic tidal 



theory it is known that the tidal radius {rt) of a globular cluster is proportional to M^^Rgc 
(Ivon Hoernerl Il957l ). where is the cluster's mass and Rgc its galactocentric distance. 
Therefore, since red globular clusters typically orbit in a stronger tidal field than blue clus- 
ters, they will be subject to increased tidal stripping and may then have smaller tidal radii. 
If the structural properties of the clusters, such as their central concentrations c, do not vary 
systematically with Rgc-, then the mean effective (or half-light) radius should also increase 
with Rgc- This statement is supported by many observa tional studies which f ind a general 



increase in cluster size w i th ga l actocentric di s tance fe.g. Ivaii den Berghl Il996l: ISpitler et al. 



2006 



Gomez fc Woodlevll2007l : iHarrisI l2009bl : IWoodlevI 120121 ). Therefore 



Larsen fc Brodie 



(120031 ) concluded that blue clusters are only larger than red clusters at all proj ected d i stances 
because they have larger three dimensional galactocentric distances. However, iHarrid (j2009bl ) 
showed that for six giant elliptical galaxies the ratio between mean red cluster size and mean 
blue cluster s i ze did not change with galactocentric distance. This w as also observed by 
Paolillo et all J201lh in NGC 1399 and lBlom. Spitler. fc Forbe^ J2012h in NGC 4365. One 
would expect that at larger galactocentric distances, proie ction effects would be n egligible 



and that red and blue clusters would have similar sizes (see ILarsen fc Brodid (120031 )). These 
findings strongly suggest that the size difference between red and blue clusters is intrinsic. 
Studies w hich support the pro jection effect hypothesis include observations of the Sombrero 



galaxy by lSpitler et al.l (120061 ) . however this is not a giant elliptical galaxy like M87. 



Explanations in the literature which propose that the size difference between red and 
blue clusters is real suggest it i s the result of red and blue clusters having different dynamical 
and stellar evolution histories. IJordanI (120041 ) suggested that stellar evolution time scales are 
dependent on metallicity such that many internal relaxation times of dynamical evolution 
would result in metal-rich clusters being small er than metal-poor clus ters. This explanation 
is supported by re cent N-body simulations by ISchulman et al.l (120121 ) of young clu sters and 
Sippel et al.l (120121 ) of clusters evolved to a Hubble time. However, a recent paper by lDowning 
(I2OI2I ) suggests instead the size difference could be attributed to black holes in blue clusters 
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heating their central regions. Finally, I Jordan et al.l ( 120051 ) and iHarria ( l2009bl ) proposed that 



metal-rich clusters may have undergone more rapid cooling and contra ction while they were 



gaseou s protoclusters, before the majority of their stars had formed. iKundu &: Whitmore 



( 119981 ) also suggested the discrepancy is due to differences in the formation of metal-poor 



and metal-rich globular clusters. 

The purpose of this study is to further explore the possibility that in M87 the observed 
size difference is due to the projection of the different spatial distributions of red and blue 
globular clusters. While this is simply one possible explanation for the observed size differ- 
ence, confirming or ruling out this possibility will contrib ute to determining whether the size 
difference is intrinsic or artificial. In a previous paper dWebb. Sills, fc HarrisI [2OI2I ) (here- 



after WSH12), we simulated a globular cluster population orbiting in the tidal field of M87 
and reproduced the observed relationship between cluster size and projected galactocentric 
distance. In our simulation, the only intrinsic difference between red clusters and blue clus- 
ters is their assumed spatial distribution. Therefore, we can compare characteristics of the 
observed and simulated red and blue cluster populations and determine if the projection of 
the red and blue cluster spatial distributions can explain the observed size difference. More 
specifically, if the size difference is truly due to the tidal truncation of differing spatial dis- 
tributions then our simulation will show that the mean red cluster size is smaller than the 
mean blue cluster size and that red clusters are smaller than blue clusters at all projected 
galactocentric distance. 



Observations 



The HST ACS/WFC Archive images we use in this study are from program GO-10543 
(PI Baltz) . The co-added composite exposures in each filter are described in detail i n 
Bird et al.l (120 lOf ). constructed through use of the APSIS software (iBlakeslee et al. [J2003h . 
While these are the same raw images used by [Madrid et al.l (120091 ) , iPeng et al.l ( 12009! ) and 
Waters et al.l ( 12009! ) . as discussed in WSH12 our images have an improved scale of 0'.'025 px~^ 
(half the native pixel size of the camera). This produced an improv ement in the eff ective 



spatial resolution of the data compared with all previous studies (see lBird et al.ll2010 
WSII12 for a more detailed description). 



and 



Details regar ding the iden tification of globular cluster candidates and the measurement 
of cluster sizes via lKingl ( 119621 ) model fitting can be found in WSII12. The final candidate list 
contained 1290 globular clusters identified with high confidence within projected distances 
Rgc < 10 kpc. To separate the population into red clusters and blue clusters, we consider the 
colour distribution of the globular cluster candidates in Figure [T] The vertical dotted line 
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in Figured] at (V-I) = 1.07 marks a local minimum in the colour distribution separating the 
blue globular cluster population (left of vertical line) and the red globular cluster population 
(right of vertical line). 

In Figure [2] we compare red and blue cluster sizes by showing the distribution of the 
average effective radius of each cluster in the V and I bands, as defined by standard [King 
(1l962l ) profiles. In general, red globular clusters are smaller than blue globular clusters. We 
find the overall mean red cluster size (2.65 pc) to be approximately 30% (0.71 pc) smaller 
than the mean blue clu ster size (3.36 pc). T he difference of 30% is slightly higher than both 
the difference found by Madrid et all (120091 ) for M87 (24%) and the mean difference of 20% 
found in the literature for all galaxy types. This difference is a result of our strict selection 
criteria when identifyin g globular cluster c andi dates (see WSH12), a s our cluster sizes are 
comparable to those of iPeng et al.l ( 120091 ) and [Madrid et al.l (120091 ). However, since the 
simulation outlined in Section 3 is tailored to fit our observational dataset, the discrepancy 
will not affect our subsequent analysis. A Kolmogorov-Smirnov (KS) test indicates that the 
Th distribution of the red and blue clusters have essentially zero (2 x 10~^^) probability of 
being drawn from the same distribution. 

A plot of median effective radius versus projected galactocentric distance is shown in 
Figure [31 Here we see that red clusters appear to be smaller than blue clusters at all radii. 



3. Simulation 



Our simulation of a globular cluster population within M87 (outlined in WSH12) allows 
for the simulated clusters to have the same characteristics as the observed globular cluster 
population. A position {R, 6, (j)), veloci ty (f^, vg, v^), mass (M), and central concentra- 
tion (the King- model c- value ( lKing||l966l )) was assigned to each cluster. Cluster mass was 



taken from the near-universal globular cluster luminosity function (e.g. iBrodie &: Strader 
20061) . with luminosities con v erted to mass by assuming a mass-to-light ratio of 2 (e.g. 
McLaughlin fc van der Marell (120051 )). The simulation allowed for clusters to have a range 
of central concentrations, with the distribution set equal to that of the Milky Way cluster 
population (Harris 1996 (2010 Edition)). 

The radial velocity dispersion of the globul ar clusters was se t equal to the line of site 
velocity dispersion of the observed M87 clusters ( ICote et al.l 120011 ). In WSH12, in order for 
the simulation to match the observed relationship between cluster size and projected galac- 
tocentric distance we allowed for outer region clusters to have increased orbital anisotropy 
such that they prefer radial orbits. This was included in the simulation presented here, and 
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Fig. 1. — Colour distribution of observed globular cluster candidates in M87. The vertical 
dotted line marks the separation between red and blue globular clusters. 
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Fig. 2. — Effective radius (r/^) distribution of observed globular cluster candidates in M87. 
The red histogram corresponds to red globular clusters and the blue histogram corresponds 
to blue globular clusters. 
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Fig. 3. — Relationship between median effective radius (r/i) and projected galactocentric 
distance (Rgc) for observed globular cluster candidates in M87. The red line corresponds 
to red globular clusters and the blue line corresponds to blue globular clusters. Each bin 
contains 25 globular cluster candidates. Error bars represent the standard error a/ \/ln) as 
given by [Harris et al.l (|2010| ) . 
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was used to calculate the tangential velocity dispersion of the simulated clusters. 



Our simulation also allowed for red and b lue clus t ers to have two separate spatial distri- 
butions. These distributions were taken from iHarrij ( l2009a[ ). who found that the projected 
radial profile of the blue and red globular cluster subpopulations could be fit with a standard 
Hubble profile relating density (ad) to projected distance (R): 



R 

aa{R) = ao/(l + — 



(1) 



The appropriate values for R^ and a are 2.0' and 1.8 for blue clusters and 1.2' and 2.1 
for red clusters. A three dimensional radial profile was then obtained by assuming that 
the globular cluster population of M87 is spherically symmetric, and integrating Equation 
[T] over a spherical volume. This results in our simulated projected radial profile matching 
Equation [T] within uncertainty. The angular distribution was also assumed to be spherically 
symmetric. 

Exactly 6000 globular clusters were simulated so the number of clusters within 10 kpc 
of M87 would be the same as the observed dataset. Clusters were simulated out to 100 kpc 
following the observed profile to ensure that larger distant clusters could be projected to 
small projected distances. 40% of the clusters were designated as "red" clusters and 60% 
were designated as "blue" clusters such t hat the ratio of blue clusters to red clusters is in 
agreement with the data in iHarrid (j2009a] ). 



Using the simulated par ameters of each cluster, we calculate the theoretical tidal radius 
of each cluster as derived by iBertin fc Varril (120081 ) : 



1/3 



(2) 



Where f2, n and v are defined as: 



= {d^a{R)/dR)RjRp 

= 3fi2 + {(f^G{.R)ldR^)R^ 
2 /ri2 



V 



4 - k'/^' 



(3) 
(4) 
(5) 



$G is the galactic potential determined from the mass profile of M87, Rp is the cluster's 
perigalactic distance, Vt is the orbital frequency of the cluster and k is the epicyclic frequency 
of the cluster at Rp. Tidal truncation is the only effect on cluster size considered in this 
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simulation, which assumes that all clusters are tidally filling. It should be noted that whether 
a cluster is red or blue has no effect on the calculation of its tidal radius. The mass profile 



of M87 was taken from iMcLaughlinl (119991 ). To best compare with our observations, tidal 



radii are c o nvert ed to effective radii by assuming each simulated cluster can be represented 



by a iKingl (119621 ) model. 



To determine if the size difference between red and blue globular clusters is due to 
projection effects we re-plot Figure [H but this time with the simulated cluster population in 
Figure HJ It should be noted that Figure S] only contains simulated clusters with projected 
distances within 10 kpc, just as the observations do. Here we see that the mean simulated 
red globular cluster size (3.31 pc) and the mean blue globular cluster size (3.36 pc) are 
essentially the same. Scatter about the mean value is simply due to the distribution of 
concentrations, velocities, and masses allowed by our model. The KS test cannot reject the 
hypothesis of the two samples being drawn from the same parent distribution with more 
than a 79% confidence level, a much stronger agreement than the observed distributions. 

While it remains true that outer region clusters experience weaker tidal truncation 
than inner region clusters, the lack of a size difference in our simulation indicates that the 
projection of larger, outer region blue clusters is not a strong enough effect to explain the 
observed size difference between the mean red and blue cluster sizes. The fact that outer 
region clusters have preferentially radial orbits will also contribute to the lack of a size 
difference in our simulation as some outer region clusters will still travel deep into the tidal 
field of M87, increasing the effect of tidal truncation. 

We further explore the influence of projection effects by plotting the relationship between 
cluster size and projected galactocentric distance for red and blue clusters in Figure [51 At 
a given projected galactocentric distance, again the median red cluster size is equal to the 
median blue cluster size. The results of our simulation are in direct disagreement with the 
observations in Figures [2] and |3l which showed that at all projected distances blue clusters 
are larger than red clusters. This statement is still true even if the model cluster population 
had an isotropic distribution of orbits. Therefore, the size difference between red and blue 
globular clusters cannot be explained by just the projection of two sub-populations that only 
differ by their spatial distributions. 



4. Summary 



While projection effects have long been one possible explanation as to why red clusters 
are smaller than blue clusters, we show that the size difference between red and blue clusters 
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Fig. 4. — Effective radius [ry^ distribution of simulated globular clusters in M87. The red 
histogram corresponds to red globular clusters and the blue histogram corresponds to blue 
globular clusters. 
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Fig. 5. — Relationship between effective radius (vh) and projected galactocentric distance 
{Rgc) for simulated globular clusters in M87. The red line corresponds to red globular clusters 
and the blue line corresponds to blue globular clusters. Each bin contains 25 simulated 



globular clusters. Error bars represent the standard error cr/ ^/{n) as given by [Harris et al. 

feoioh . 
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cannot be explained by projection effects alone, and that another mechanism is playing a 
more dominant role. In fact, this work suggests that red clusters must either under-fill their 
tidal radius or have different orbital parameters than the blue population in order to be 
sufficiently smaller than blue clusters. Specifically, if red clusters had preferentially radial 
orbits compared to blue clusters, it would result in red clusters having smaller tidal radii than 
blue clusters at all galactocentric distances. We suggest that the metallicity of a globular 
cluster must play an important role in determining exactly how large or small a globular 
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Applying our model to other galaxies with different galactic potentials and globular 
cluster systems will confirm whether this conclusion applies to all galaxies, all early-type 
galaxies, or just M87. Future work will explore the assumptions made in our simulation, 
including the relationship between tidal radii and effective radii and additional differences 
between the red and blue clusters in our model, such as the possibility of the sub-populations 
having different orbital anisotropy profiles. It is also worth exploring what influence globular 
clusters that do not fill their tidal radius have on our comparison, as some studies find 
cluster sizes increase at a much shallower rate at large galactocentric distances (or not at 



all) compared to o ur model (e.g. iGomez fc Woodleyl 120071 ; ISpitler et al.ll2006t lHarrisll2009b 



Harris et al.l 120101). Finally , simulations of clusters with a range of metallicities, similar to 



those of ISippel et al.l (120121 ) which find the size difference can be attributed to the effect that 
metallicity has on stellar evolution and dynamical effects, will also shed a significant amount 
of light on this issue. 
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